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Abstract 

Recently, a general analysis has been given of the stability with respect to axisymmetric pertur¬ 
bations of stationary-axisymmetric black holes and black branes in vacuum general relativity in 
arbitrary dimensions. It was shown that positivity of canonical energy on an appropriate space of 
perturbations is necessary and sufficient for stability. However, the notions of both “stability” and 
“instability” in this result are significantly weaker than one would like to obtain. In particular, if 
there exists a perturbation with negative canonical energy, “instability” has been shown to occur 
only in the sense that this perturbation cannot asymptotically approach a stationary perturbation 
at late times. In this paper, we prove that if a perturbation of the form £thg —with 6g a solution to 
the linearized Einstein equation—has negative canonical energy, then that perturbation must, in 
fact, grow exponentially in time. The key idea is to make use of the t- or (t-^)-reflection isometry, 
i, of the background spacetime and decompose the initial data for perturbations into their odd and 
even parts under i. We then write the canonical energy as S' = JC + '^, where JC and ‘Hi, respec¬ 
tively, denote the canonical energy of the odd part (“kinetic energy”) and even part (“potential 
energy”). One of the main results of this paper is the proof that JC is positive definite for any black 
hole background. We use JC to construct a Hilbert space Jif on which time evolution is given in 
terms of a self-adjoint operator A, whose spectrum includes negative values if and only if fails 
to be positive. Negative spectrum of A implies exponential growth of the perturbations in Jif that 
have nontrivial projection into the negative spectral subspace. This includes all perturbations of 
the form £tdg with negative canonical energy. A “Rayleigh-Ritz” type of variational principle is 
derived, which can be used to obtain lower bounds on the rate of exponential growth. 
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1. INTRODUCTION 

Determining the stability of exact solutions to Einstein equations is a long standing 
problem. In order for a solution to be physically relevant, it is essential that sufficiently small 
perturbations not drive one away from that solution. However, the issue of full nonlinear 
stability is extremely difficult to analyze and has been settled only for Minkowski spacetime 
[1]. As a hrst step, it is important to analyze the stability of solutions to linear perturbations, 
which solve the linearized Einstein equations (or in general, other helds that satisfy a linear 
dynamical equation). 

A class of solutions of considerable interest are black hole and black brane spacetimes. 
All of the results of this paper will apply to the linear stability with respect to axisymmet- 
ric perturbations of stationary and axisymmetric black holes and black branes in vacuum 
general relativity in arbitrary dimensions. In order to keep our analysis and notation less 
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cumbersome, we shall restrict consideration in this paper to the case of black holes, but the 
generalization of our results to the case of black branes is straightforward. 

In 4-spacetime-dimensions, the Schwarzschild and Kerr solutions have long been known 
to be the unique stationary black hole solutions to the the vacuum Einstein equation. These 
solutions are believed to physically describe the asymptotic hnal states of gravitational 
collapse. In order for this to be the case, it is essential that these solutions be stable. In 
higher dimensions, there exist many other types of black hole solutions, which are of interest 
for various theoretical reasons. It is of considerable interest to analyze the stability of these 
solutions as well. 

To establish the linear stability of a solution, one must show that all suitably regular 
(e.g., smooth and satisfying appropriate asymptotic conditions) initial data for the linearized 
equations give rise to solutions that remain uniformly bounded in time. In the case of a 
linearly stable, stationary black hole solution, one would like to establish an even stronger 
result, namely that all suitably regular linearized solutions decay at asymptotically late 
times to a stationary solution. On the other hand, a considerably weaker notion of stability 
that is much easier to analyze is mode stability, i.e., the nonexistence of suitably regular 
solutions that grow exponentially in time. Clearly, mode instability implies instability in 
any other reasonable sense, but mode stability does not directly imply uniform boundedness 
or decay properties of linearized solutions. 

As a simpler, model problem, it is useful to study scalar held perturbations of a black 
hole. For the case of a Schwarzschild black hole, it is easy to prove mode stability. Hilbert 
space methods were used in [2] to establish stability of Schwarzschild in the sense of uniform 
boundedness of solutions with a given spherical harmonic angular dependence. However, in 
addition to the restriction to a given spherical harmonic angular dependence, this approach 
required the imposition of the unwanted restriction of the vanishing of the perturbation on 
the bifurcation 2-sphere. This latter restriction arose from the requirement that the pertur¬ 
bation lie in the Hilbert space used in the analysis and, as pointed out in [3], it is equivalent 
to the requirement that the perturbation be expressible as Tj of another perturbation, where 
£t denotes the Lie derivative with respect to the timelike Killing held, t^, of the background 
spacetime. The restriction to a given spherical harmonic angular dependence can be straight¬ 
forwardly removed and, indeed, the method can be generalized to analyze the stability of 
any static spacetime to scalar perturbations [4]. However, it required a “trick” for Kay and 
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Wald [5] to remove the restriction of the vanishing of the perturbation on the bifurcation 2- 
sphere, thereby establishing uniform boundedness of scalar perturbations of Schwarzschild. 
Nevertheless, the methods used in [2, 4, 5] are inadequate for establishing decay of per¬ 
turbations. More recently, improved methods, described in [3], have established decay of 
a scalar held in Schwarzschild and, indeed, similar results have been proven to hold in the 
much more difficult case of a Kerr black hole (including non-axisymmetric perturbations) 
[6-9]. However, these methods rely on very detailed properties of the Schwarzschild and 
Kerr metrics and cannot straightforwardly be generalized to arbitrary higher dimensional 
black holes. 

For the case of gravitational perturbations, mode stability of 4-dimensional Schwarzschild 
spacetime follows immediately from the form of the decoupled equations originally obtained 
by Regge and Wheeler [10] and Zerilli [11]. These results have been extended to prove mode 
stability of higher-dimensional Schwarzschild spacetimes by Ishibashi and Kodama [12]. The 
methods of [2, 5] could then be applied to show uniform boundedness of solutions with a 
given spherical harmonic angular dependence. However, on account of the gauge choices 
made, it is not straightforward to remove the restriction to a given spherical harmonic 
angular dependence. The decay results of [3] have not yet been generalized to gravitational 
perturbations, even for the case of Schwarzschild. 

Recently, a new approach to investigating the linear stability of arbitrary static or sta¬ 
tionary and axisymmetric black holes and black branes to axisymmetric gravitational per¬ 
turbations was given in [13]. A conserved, symmetric, gauge-invariant quadratic form on 
initial data for perturbations, called the canonical energy, was introduced and shown to be 
given by the formula 

^ = 5‘^M - - —6'^A ( 1 . 1 ) 

A 

where M denotes the ADM mass, Ja denotes the ADM angular momenta (with the sum 
over A being taken over independent “rotational planes”), and A denotes the surface area 
of the horizon. For perturbations with 6M = 5 J\ = 5Pi = SA = 0 (where Pi denotes the 
ADM linear momentum), it was shown that S is non-degenerate on the space of linearized 
solutions modulo gauge and modulo perturbations to other stationary and axisymmetric 
black holes. Consequently, either S is positive definite on this space or it can take negative 
values. Positive dehniteness of S immediately establishes mode stability, since the existence 
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of an exponentially growing solution is incompatible with the presence of a positive dehnite 
conserved norm. On the other hand, if S’ can take negative values, then a further argument 
establishes instability as follows: It was shown in [13] that for axisymmetric perturbations, 
the flux of S through null inhnity and through the black hole horizon must be positive. 
Thus, if the perturbation were to suitably approach a perturbation to another stationary 
black hole at asymptotically late times (with the limit suitably taken along the orbits of 
the timelike Killing held of the background), then S' < S, where S' denotes the canonical 
energy of this perturbation to another stationary black hole. Consequently, if < 0, then 
S' < 0. Furthermore, for the limiting perturbation, we continue to have dM = (5Ja = 0, 
since huxes of mass and angular momenta are quadratic in the perturbation. But S' < ^ 
then contradicts the fact that the canonical energy is degenerate on perturbations to other 
stationary black holes with 6M = 5Ja = 0, thus establishing “instability” in the sense that 
the perturbation cannot asymptotically approach a perturbation to another stationary black 
hole. 

However, the analysis of [13] leaves many questions unanswered with regard to the sta¬ 
bility and instability of black holes to axisymmetric gravitational perturbations. In the 
hrst place, one would like to know—for any given black hole—whether the expression for 
S is positive. An explicit integral expression for S in terms of initial data qab) for the 
linearized perturbation is given in Eq.86 of [13], where Qab denotes the perturbation of the 
spatial metric and Pab denotes the perturbation of the ADM momentum variable. Thus, 
the issue is simply to determine, for a given background, whether this expression is positive 
for all perturbations. However, the expression for S is quite complicated, and {pab, qab) are 
not “free” but must satisfy constraints and boundary conditions. Furthermore, it would be 
difficult to show that the integral expression for S is positive without writing it in a form 
where the integrand is positive; however, although S is gauge invariant, the integrand is 
not, so this would require an appropriate choice of gauge and it is not obvious what gauge 
conditions will work. (The gauge conditions we will introduce later in this paper do not 
seem to work for this purpose.) For the case of a thermodynamically unstable black brane, 
there is a relatively obvious candidate perturbation that makes S negative, and failure of 
the positivity of S in this case was proven in [13]. However, establishing the positivity of 
S in cases where it should be positive seems quite difficult. Indeed, even for the case of 
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Schwarzschild, we have been able to prove positivity of S’ only by an indirect^ argument— 
see remark 3 following Thm. 1 of Sec. 4 below. We will not consider further in this paper 
the issue of determining the positivity of S in specihc spacetimes. 

As already mentioned above, positivity of S immediately implies mode stability, but if 
S is positive, one would expect signihcantly stronger results to hold, namely, the decay of 
perturbations. It is possible that the methods introduced in this paper may be adequate 
to prove uniform boundedness of perturbations when S is positive, i.e., the methods of this 
paper may be adequate for generalizing the results of Kay and Wald [5] to axisymmetric 
gravitational perturbations of general stationary-axisymmetric black holes and black branes 
for which S is positive. However, considerable further analysis would be required to show 
this^ and we shall not attempt to carry out this analysis here. Our methods would not, in 
any case, be adequate for proving decay results. To prove decay results, one would like to 
show that a suitably modihed version of S is “coercive.” It would seem that a necessary 
hrst step toward showing this would be to have a good technique to show that S is positive 
(in cases where it is positive), but, as discussed in the previous paragraph, this is currently 
lacking. We shall not consider further in this paper this issue of obtaining strengthened 
stability results when S is positive. 

As explained above, if S fails to be positive on a black hole or black brane spacetime, 
then it was shown in [13] that one has instability in the sense that there exist perturbations 
that cannot approach a stationary perturbation at late times. However, if S fails to be 
positive, one might expect a much stronger result to hold, namely, the existence of initially 
well behaved perturbations that grow exponentially with time. The main purpose of this 
paper is to prove that this is indeed the case. Specihcally, we will prove that if a perturbation 
of the form £t^g has negative canonical energy, then that perturbation must, in fact, grow 
exponentially in time^ (in all gauges). We now outline the key ideas used in the proof. 

The hrst key step—undertaken in Sec. 3—is to completely hx the gauge of our pertur- 
^ In the case of Rindler spacetime, we have succeeded in proving positivity of S directly from the formula 
for S. 

^ The main things that would have to be shown are that (i) the operator A of Eq. 5.16 and its powers provide 
norms that, together with our gauge conditions, are equivalent to Sobolev norms and (ii) a version of the 
Kay and Wald “trick” can be used to eliminate the restriction to perturbations of the form of £t applied 
to another perturbation. 

^ In the last sentence of [14], the authors raise the question of whether violation of the local Penrose 
inequality implies the existence of perturbations 6hat grow exponentially with time. Since it was shown 
in [13] that violation of the Penrose inequality is equivalent to the failure of S to be positive, our results, 
in essence, answer that question in the affirmative. 



bations so as to obtain a unique time evolution. We consider a maximal slice S of the 
background spacetime"^ and work in the space, of linearized initial data 6X = (pab, gab) 
on S that lie in the intersection of the weighted Sobolev spaces dehned by Eq. 3.3 be¬ 
low. Following the strategy of [13] we impose the linearized constraints, the conditions 
6M = 5J\ = SPi = 0, and certain gauge conditions at the horizon by means of a projection 
map flc. We then fix the gauge completely by applying a similar projection map 11^ that 
commutes with lie. The needed properties of lie and 11^ are proven in Appendix A. The 
space contains a gauge representative of any SX G that satisfies the linearized 

constraints and the conditions SM = 5J\ = 0. The space := 'n.gYlc[^'^] contains a 
unique gauge representative of any such 5X, and it therefore provides a suitable space to 
study dynamics. 

The next key idea is to make use of the existence of a suitable reflection isometry, i, of 
the background spacetime. In the case of a static black hole, the maximal slice S must be 
orthogonal to the static Killing held P. The desired isometry, i, can be constructed by “f- 
rehection” about S, i.e., i is the diffeomorphism that takes a point p lying at proper time r 
along a normal geodesic 7 starting from g G S to the point p' lying at proper time —r along 
7 . For the case of a stationary-axisymmetric black hole, it has recently been proven [16] that 
if the stationary-axisymmetric isometries act trivially^, then there exists a “(f-^j-rehection” 
isometry, i, dehned as follows®: Let (j)'^,..., denote the collection^ of commuting axial 
Killing helds. Then, as shown in [16] the axial Killing helds are tangent to S and are surface 
orthogonal within S, so there exists an isometry, is : F —)■ S, that rehects about a surface, 
S, in S. Now, extend is to a spacetime difeomorphism i by mapping a point p lying at 
proper time r along a normal geodesic starting from g G S to the point p' lying at proper 
time —r along the normal geodesic starting at is{q)- Then, as shown in [16], i preserves the 
initial data on S and thus dehnes an isometry. We restrict consideration in the remainder 
of this paper to black holes that possess a t- or (t-^j-rehection isometry i. 

Existence of S is guaranteed by the results of [15]. 

^ By “act trivially,” we mean that if we remove from the spacetime manifold the points at which the 

stationary-axisymmetric Killing fields are linearly dependent, then the resulting manifold acquires the 

structure of a trivial principal fiber bundle with respect to the stationary-axisymmetric symmetries. This 

rules out behavior of axial Killing fields similar to that occurring in the Sorkin monopole [17] [18]. 

® This result is well known in 4-spacetime dimensions (see, e.g., [19]), but, as explained in [16], the proof 

in 4-dimensions does not generalize to higher dimensions, so the additional assumption of a trivial action 

is needed in spacetime dimension greater than 4. 7 

^ We may, if we wish, delete from this collection any axial Killing fields with vanishing horizon rotation 
provided that the action of the remaining Killing fields is also trivial [16]. In particular, for a static black 
hole, we may delete all of the axial Killing fields and take i to be the “t-reflection” isometry, as assumed 



The next key step—undertaken in Sec. 4—is to decompose linearized initial data 5X = 
{PabjQab) ^ luto its odd and even parts under the action of i, i.e., we decompose the 
space, of allowed initial data for perturbations as W^e denote 

elements of by “P” and we denote elements of by “Q,” so, by construction, we 

have i*P = —P and i*Q = Q, where i* denotes the action induced by i. In the static case, 
we have P = (pafe,0) and Q = {0,qab)- In the stationary-axisymmetric case, P is a linear 
combination of the “polar” part of Pab and the “axial” part of Qab, whereas Q is a linear 
combination of the “polar” part of qab and the “axial” part of Pab (see Eq. 4.2 below). The 
linearized constraint equations are automatically invariant under i. Our gauge and boundary 
conditions are also chosen to be invariant under i, so P and Q decouple and we may treat 
them as independent perturbations. The decomposition y°^ = y^^ © allows us to 
break up the canonical energy S’ into a sum of “kinetic” and “potential” energy, 

S = jy + ^ ( 1 . 2 ) 

where jy is the restriction of S to and ^ is the restriction of S to no “(P- 

(5)-cross-terms” can arise on account of the reflection symmetry. Note that XS dehnes a 
symmetric quadratic form on y^^ and ^ dehnes a symmetric quadratic form on 

The next key result—obtained at the end of Sec. 4—is the proof that jy is positive 
dehnite on i.e., 

jr(P,P)>0, (1.3) 

with equality holding only if P = 0. In particular, this shows that any failure of positivity 
of S must arise from the “potential energy,” ^. 

Time evolution is considered in Sec. 5. The ADM evolution equations mix the odd and 
even parts of the perturbation. Specihcally, for axisymmetric perturbations, the evolution 
equations take the general form 


Q = /CP 


(1.4) 


P = -UQ 


(1.5) 


for some operators /C : —)■ lA : y^di where the overdot denotes the 

Lie derivative, £t, with respect to P. Furthermore, the operators /C and U appearing in the 
evolution equations are related to the above kinetic and potential energy quadratic forms 



by 


JY{P,P) = Qj,{P,}CP) 

(1.6) 

^{Q,Q) = -nj,{Q,i(Q) 

(1,7) 


where fis denotes the symplectic form. These equations express the fact that the canonical 
energy S is actually a Hamiltonian for the linearized system. 

As we show in Sec. 5, the positivity of allows us to define a “/C“^ Hilbert space,” 
as follows. Consider Q of the form Q = /CP for some P; Q will be of this form if and only 
if the initial data P' = 0, = Q corresponds to a linearized solution that can be expressed 

as £t of another linearized solution. Dehne the inner product oi Q = /CP and Q = /CP by 

(Q,g) :=^(P,P) (1.8) 

and dehne to be the Hilbert space completion of this inner product space. 

Combining the two time evolution equations Eq. 1.4 and Eq. 1.5, we obtain the following 
second order evolution equation involving Q alone: 

Q = -AQ (1.9) 

where 

A := KIA (1.10) 


is a symmetric operator on the Hilbert space Consequently, by passing to a self-adjoint 
extension of A, we can solve Eq. 1.9 by spectral methods. For initial data in fl we 
show that the Hilbert space solution must coincide with the PDE solution. It then follows 
that if there exists initial data of the form Q = /CP for P G for which {Q^AQ) ^ = 
'Y{Q,Q) < 0, then the ^-norm of Q must grow exponentially with time®. The gauge 
invariant quantities JY and for this perturbation must also blow up exponentially, thus 
showing that the exponential blow up is not a gauge artifact. 

The above instability result—including a quantitative statement about the rate of expo¬ 
nential growth—can be formulated as a “Rayleigh-Ritz variational principle” as follows. Let 
P G Y^^ be any reflection-odd smooth initial data satisfying the constraints and such that 
SJj, = SPi = 0. Write 


2 _ ^(/CP,/CP) 
“ jr(p, p) 


( 1 . 11 ) 


The earliest reference that we are aware of for the basic argument that an equation of the form of Eq. 1.9 


will have exponentially growing solutions if A has negative spectrum is [20]. 
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Then if < 0, the solution generated by {P,Q = 0) will grow at least as fast as exp(|a;|t) 
(in all gauges). 


The remainder of this paper is devoted to fleshing out these arguments. Sec. 2 discusses 
the properties of the background black hole spacetime and introduces the ADM formal¬ 
ism that we use throughout. In Sec. 3 we introduce the spaces of interest for linearized 
perturbations and define the projection maps lie and 11^, which completely gauge-fix the 
perturbations. Sec. 4 introduces the splitting of the canonical energy into kinetic and poten¬ 
tial energies and ends with the theorem showing that the kinetic energy is positive-dehnite. 
In Sec. 5 we formulate the ADM equations as dynamical evolution equations on the Hilbert 
space Jif. In Sec. 6 we solve these equations by spectral methods and prove that if the 
potential energy can be made negative on elements of then there exist linear pertur¬ 
bations which have exponential growth in time. We also obtain our variational principle 
formulation. Appendix A provides the details of the construction of the projection operators 
He and H^. Appendix B provides a parallel analysis for the Klein-Gordon scalar held and 
the electromagnetic held on a hxed black hole background. 


In this paper, lower case Greek indices will be used to denote tensors on spacetime, e.g., P 
denotes the timelike Killing held of the background black hole. Lower case Latin indices will 
be used to denote tensors on the initial data surface S, e.g., qab denotes the perturbed metric 
on S. Upper case Latin indices will be used to denote tensors on the bifurcation surface, B, 
of the black hole, e.g., denotes a vector held on B. We will use the index “r” to denote 
projections normal to B, e.g., if Va is a one-form held on S, the normal component of its 
restriction to B will be denoted Vr = r°‘Va, where r“ is the unit normal to B. The spacetime 
derivative operator of the background black hole will be denoted as the background 
derivative operator on S will be denoted as Da'-, the background derivative operator on B 
will be denoted as If we consider the restriction to S of a spacetime vector held it 
will often be useful to represent as the pair (^,^“) where ^ := (with the unit 

normal to S) is the component of normal to S while is the projection of tangent to 
S. We will use the notation ^ when we wish to view in this way. 
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2. BACKGROUND SPACETIME 


We consider a (d + l)-spacetime-dimensional, asymptotically flat, static or stationary- 
axisymmetric black hole spacetime {M,g) shown in Fig. 1, with a bifurcate Killing horizon 
H := U H~, and bifurcation surface B := fl H~. Let denote the time translation 
Killing held, i.e., the Killing held that is timelike near infinity. Let S be an asymptotically 
flat Cauchy surface for one of the exterior wedges, smoothly terminating at B. Below, we 
will choose S to be a maximal slice but we need not make this choice now. We also assume 
that the bifurcation surface B is compact, but we do not assume any further restrictions 
on its topology. Let St denote the foliation obtained by applying time translations to S. 
Let denote the future-directed unit normal to S. We decompose into its normal and 
tangential parts relative to S, referred to as the lapse, N = and shift, iV“, on S. 



FIG. 1. Carter-Penrose diagram of the black hole spacetime {M,g). 


Let hab denote the induced metric on S, with determinant with respect to some fixed 
volume form denoted by h, and let Kab denote the extrinsic curvature of S. The phase space 
of general relativity is the set of initial data on S denoted by X = hab) where 7r“^ : = 

y/h(K^b _ ^ f^ab^ 

is the canonical momentum-density conjugate to hab with K := Kabh^"^. 
In order to correspond to a solution to Einstein’s equation, the initial data must satisfy the 
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constraint equations = 0 , with 


(£ =-i? + i (7r“Vafe - ^TT^) (2.1a) 

r ^ -20, ( J) (2,1b) 

where tt := 7 rafe/i“^, Da denotes the covariant derivative compatible with hab, and Rab is the 
Ricci curvature of hab- 

On the bifurcation surface R, we introduce a unit normal vector r“ (pointing into S) and, 
without loss of generality, extend r“ in a neighbourhood of B such that it is a geodesic, so 
that, in particular, {r‘^Darb)\B = 0. The projector onto B is Sab '■= hab — ^'a^h- As noted at 
the end of the introduction, we will use capital Latin letters A, B,C,... to denote tensors 
on B and an index r to denote projections normal to B. The induced metric on B will be 
denoted as sab- We denote the metric-compatible covariant derivative operator on B by 
and we write dr '■= r'^Da- The induced volume element on B will be denoted as £Ai...Ad_i- 
Note that N\b = 0 and {r^DaN)\B = k, where k is the surface gravity of the black hole 
(assuming that the axial Killing helds are tangent to S). We also have ("dijls = 0, where 

:=s^^{Kab±Darb) , ( 2 . 2 ) 


i.e., ('d+)|B is the expansion of the outgoing null geodesics at B and ('d_)|B is the expansion 
of the ingoing null geodesics at B. 

The asymptotic flatness conditions on our stationary black hole are that there exist 
coordinates (xi,..., Xd) on S such that 


hab~5a6 + 0(l/p"-2) ; 7r“' ~ 0(l/p''-') 

TV ~ 1 + 0{l/p'^-^) ; iV“ ~ 0(l/p'^-2) 


(2.3) 


where p = (xf -f ... - 1 - near inhnity. In addition, derivatives of the above quantities 

are required to fall off faster by an additional factor of 1/p^. The asymptotic conditions on 
the lapse and shift ensure, in particular, that D goes to an asymptotic time-translation at 
inhnity. 

The ADM time evolution equations are (see, e.g., Sec.E.2 of [19] for d = 3 and Sec.VI .6 
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of [21] for general d): 

1 -ab ]\T f Tjab 1 TDi.ab\ , N i^ab f^cd^ 1 ^ 2 \ 2N f^ac^ b 1 ——ab 

^TT =-N[R -^Rh ) + ^h [tT Tied- ^TT j - — (tT TTe - ^TTTT 

+ D‘^D'^N - h^^AN + De ^ 7r^‘^^DeN^'> 

2N 

hab = {^ab - ^TT/lafe) + ^D^aNb) (2.4b) 

where the overdot denotes £t and A := D°‘Da is the Laplacian on S. Since we are considering 
stationary black hole spacetimes, the left side of Eq. 2.4 vanishes. 

We can signihcantly further simplify the right side of Eq. 2.4 by choosing S to be a 
maximal slice, whose existence® was proven in [15]. As shown in [16], the axial Killing 
helds with A = 1,... ,p, are then tangent to S, so we also may denote them as 4)%. 

We will assume the existence of a (t-(;/))-reflection isometry about S, which, as discussed 
in the introduction, has been proven to exist (see [16]) if the action of the stationary and 
axisymmetric isometries is trivial. At a point of S, the (t-(;/))-reflection isometry reverses the 
directions of and but leaves the subspace orthogonal to and invariant. Since the 
normal, to S must reverse sign under the isometry, it must be expressed as a (variable) 
linear combination of and the axial Killing helds (p^,... ,(p^ 

u^{x) = a{x)[t^— ''^‘^N^{x)(p4P\ ■ (2.5) 

A 

Thus, the shift vector takes the form 

AT" = ^ N^(Pl . (2.6) 

A 

The extrinsic curvature of S is purely “axial” (i.e., odd under ^-rehection), so takes the 
form 

= (2-7) 

A 

with 7r^00 = 0 for all A, 0. 

The axial Killing helds comprise a vector space V, and it is useful to think of the A-index 
in as an abstract index of V rather than a labeling index running from 1 to p. At each 
a; G S where the Killing helds are linearly independent, we can then dehne a positive dehnite 
inverse metric <hAe(a:) on V by 

^Ae{x) ■= hab{x)(pl(pQ ( 2 . 8 ) 

® The proof of [15] was given in the case of 4 spacetime dimensions, but it generalizes straightforwardly to 
arbitrary dimensions. 


(2.4a) 
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We will use <I>a 0 and its inverse, to lower and raise V-indices. Note, however, that 
Da^Ae 7 ^ 0 so, while 0^ satishes Killing’s equation, 0“^ does not. The (d — p)-dimensional 
surface orthogonality of 0^ within S together with Killing’s equation implies that 

Da4>bA = —*h®“00[aT>fe]‘hAE: • (2.9) 

With the above choice of S, the ADM evolution equations Eq. 2.4 reduce to: 

NRah = DaDbN - 2N {iT^TTbA - 7r^7r^®0aA0fe0) (2.10) 

and 

DaN^ = -2N7T^ 

In addition, the constraint equations Eq. 2.1 become: 

R = 

and 

DaK = 0 (2.13) 

These relations simplify considerably in the static case, where 7 r“^ = 0 and = 0. The 
evolution equations then reduce to 


( 2 . 11 ) 

( 2 . 12 ) 


NR,b = DaDbN (2.14) 

and the constraint equations reduce to 

i? = 0. (2.15) 

3. LINEAR PERTURBATIONS: CONSTRAINTS, BOUNDARY CONDITIONS, 
AND GAUGE CONDITIONS 

Let X(A) = ( 7 r“'’(A), hab(A)) be a one-parameter family of initial data that is jointly 
smooth in A and point on S, with X(0) corresponding to initial data for a stationary 
black hole, as discussed in the previous section. Linearized perturbations off of W(0) are 
characterized by 

6X = {pab,qab) (3.1) 
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where 


yhp«'>;=57r“' = -^7r“'(A)U=o ; Qab ■= SKb = 4^Kb{\)\x=o (3.2) 

d\ d\ 

The main task of this section is to dehne the space of perturbations, that we will 
work with. We wish to dehne a space of perturbations that include all perturbations that (i) 
are smooth and satisfy appropriate asymptotic fall-off properties at inhnity, (ii) satisfy the 
constraints, (iii) satisfy the horizon gauge conditions needed to dehne the canonical energy, 
and (iv) have vanishing perturbed ADM mass, linear momentum, and angular momentum 
as needed for the stability analysis. Furthermore, we want the perturbations to satisfy the 
property that (v) they are completely “gauge hxed.” 

The smoothness and asymptotic fall-oh properties will be enforced by requiring the initial 
data 5X on S to he in suitable weighted Sobolev spaces. Let p be a positive function that 
goes to 1 in a neighborhood of B and approaches {x\ -|- ... -1- near inhnity. Let 

denote the closure of data in that vanish in a neighborhood of inhnity in the norm 

k „ 

II '5V \\w^ '= E / .r>.„pfc)(c“',,, D-p*') + (£>.....Da.fcjir'-... c““P)] 

n=0 

(3.3) 

Here and below, the integral over S is taken with respect the volume form induced by the 
background metric hab- The space of interest for our analysis is 

^^:=nkWl;. (3.4) 

The family of IT^-norms gives the natural structure of a Frechet space. The hniteness 
of all the above Sobolev norms together with the standard Sobolev estimates, implies that all 
elements of are smooth, with Pab,<lab = 0 ( 1 /^*^/^) as p —)■ cx) and all spatial derivatives 
falling faster by corresponding powers of p. Note that the fall-oh conditions for data in 
are weaker than normally assumed (see Eq. 2.3 above), except for d = 3,4 where our weighted 
Sobolev space conditions impose a faster than normal fall-oh in Qab and thereby exclude 
perturbations that change the ADM mass. Since we will be interested only in perturbations 
for which SM = 0, this will not impose any unwanted restrictions when d = 3,4. Note also 
that we impose the same fall-oh conditions on pab and qab rather than requiring faster fall-oh 
on Pab- It is important for our constructions below that we treat pab and Qab on an “equal 
footing.” The imposition of the same fall-oh rates on pab and Qab would cause difficulties 
if we wished to consider time evolution to “boosted slices.” However, when we consider 
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time evolution of perturbations in Sec. 5, we will evolve only to “time translated” Cancliy 
snrfaces. 

Althongli our interest is in perturbations 6X G it is very convenient to perform 

constructions in the larger L^-space 


:= < 3 <^^g7 (3.5) 

with inner product 

(^6X, SX^ = ^ (rtafe + , (3.6) 

since, as we shall see below, the constraints and gange conditions can be expressed in terms 
of orthogonal projection maps on Note that the symplectic form 

fis (X; <5X) ;= ^ - q^bP^^) (3.7) 

is represented on J^gj- by the bounded linear map S given by 

S {Pab, Qab) = {Qab, -Pab) (3.8) 

Note further that S* = —S and = —1, so S : iX’g^ —)■ ^gr is an orthogonal map. 

Of course, we are interested only in the elements 5X G that satisfy the linearized 
constraints. We may view the operator obtained by linearizing the constraints Eq. 2.1 off of 
the stationary black hole backgronnd as a map, £, taking smooth initial data on S to a pair 


^ = K.r) 


(3.9) 


consisting of a smooth scalar held, and smooth vector held, on S: 

2 {1 ^ UCLh\ ^ I 2 {„ac„ b 1 ^^ab\ ^ \ 

^ d-l^ ^ ) Pab + d-1^^ ) 

- i - ^TT^] q - D^D\ab + Ag + R^\ab 

where q := qabh^^- In this notation, the linearized constraints are 



(3.10) 


c ;= C5X = 0. 


(3.11) 
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The formal L^-adjoint of £ is the map 

^ (j^ab dZi'^hab) + 2Z)(a'Cb) 


a 




\ 


{ (fi-“ - i - Jry) + f 

- D‘D''(l + - £>„ (f ^ 


(3.12) 


where we have used the background constraints Eq. 2.1 in computing the adjoint. Note 
that S*C*^ corresponds precisely to the inhnitesimal gauge transformation generated by the 
vector held Note also that since all gange transformations are solutions to the 

linearized constraints, the equation CS*C* = 0 holds as an identity. 

As explained in [13], in addition to the constraints, we also must impose some restrictions 
on the pertnrbations and some gange conditions at B in order that the canonical energy 
have snitable gange invariance and non-degeneracy properties. The restrictions we impose 
are 5M = 6J\ = 0, where J\ denotes the angnlar momenta conjugate to the axial Killing 
helds (j)^] we also impose SPi = 0 where Pi denotes the linear momenta, but this is not a 
physical restriction since this condition can be achieved via an asymptotic Lorentz boost. 
The gange conditions at B that were imposed in [13] were = 0 and (5'd+ = 0, where 


fe = (3.13) 

is the perturbed area element^^ of B and is the pertnrbed outgoing expansion of B. 
However, for onr pnrposes, it is essential that onr gauge conditions at B respect the (t- 
0)-reflection isometry i. For this reason, we impose the additional gauge condition that 
the ingoing expansion of B also vanish, d'd- = 0. That both = 0 and d'd- = 0 can 
be achieved withont imposing any physical restrictions on the pertnrbation can be seen as 
follows: It was proven in [13] that the gauge condition (5-d+ = 0 always can be imposed. 
However, the condition = 0 does not nniqnely determine a 2-snrface; rather it holds 
on all cross-sections of an ontgoing nnll hypersnrface in the perturbed spacetime. The 
condition = 0 similarly holds on all cross-sections of an ingoing nnll hypersnrface. The 
intersection of these hypersnrfaces dehnes a nnique snrface on which both = 0 and 


Since 6M = (5Ja = 0, the first law of black hole mechanics implies 5A = 0, in which case fc = 0 can be 
imposed by a gauge choice. 
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Si}- = 0. An infinitesimal differmorphism that moves this surface to B achieves the desired 
gauge condition. 

Thus, the conditions we impose on perturbations are 


SM = SJa = SPi = 0 (3.14) 

5e = 5d+= = 0. (3.15) 


As already explained above, only the conditions SM = SJa = 0 are physical restrictions on 
the perturbations; the remaining conditions SPi = Se = = Si)- = 0 can be achieved 

by a choice of gauge. The relations Eq. 3.15 can be written more explicitly in terms of the 
perturbed initial data as 

0 = fe = (3.16a) 

0 = S'Sodd = -Prr (3.16b) 

0 = (5deven = | (3.16 c) 


where Si)± = Si)odd ± SHeven- Here we have used the subscripts “odd” and “even” to denote 
parts of the perturbed expansion which are odd and even under the t-(f) reflection isometry, 
as will be discussed further in Sec. 4. 

Following the strategy of [13], we impose the linearized constraints Eq. 3.11 and the addi¬ 
tional conditions Eq. 3.14 and Eq. 3.15 by the following procedure. Let #0 be the subspace 
of comprised by smooth gauge transformations that become an asymptotic translation 
or a rotation with respect to an axial Killing held^^ at inhnity and whose projection onto B 
vanishes 


Wc := G I = 0 and goes to an asymptotic translation plus 

rotation with respect to any axial Killing held at inhnity} 


(3.17) 


where, as above, ^ denotes the pair consisting of a smooth scalar held, and smooth vector 
held, on S. For any ^ as in Eq. 3.17 and any SX G a direct computation yields [13] 

{C*^, SX) = (^, CSX) + [ [-2dr^ Se + ^ {S^+ + S^-) + (5^9+ - (5d_)] + SH^ (3.18) 

Jb 


Note that all axial Killing fields of the background spacetime are included in the definition of Wc- In 
particular, for perturbations of Schwarzschild, may approach an arbitrary asymptotic rotation at 
infinity. 
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where SH^ is the boundary term at inhnity that represents the perturbed ADM conserved 
quantity corresponding to the asymptotic symmetry Thus, SX G satishes the 
constraints Eq. 3.11 and the conditions Eq. 3.14 and Eq. 3.15 if and only if it is L^-orthogonal 
to C*C, for all ^ as in Eq. 3.17. 

Let Xc denote the closed subspace of ^gr that is symplectically-orthogonal to Wc, 

% := ^ I ^ Q ^ ^ 

Then elements of % weakly satisfy the constraints Eq. 3.11 and the conditions Eq. 3.14 and 
Eq. 3.15. The perturbations of interest for our analysis are those that he in the subspace 

;= % n . (3.20) 

Let Lie : ^gr —t ^gr be the orthogonal projection operator onto %. By slight modiheations 
of the proof of Lemma 3 of [13], we show in Appendix A that—by virtue of the fact that 11^ 
is constructed by solving an elliptic system—we have 

. (3.21) 

Furthermore, we show that Lie : t is continuous in the natural (Frechet) topology 

of Note also that since is dense in in the L^-topology and Lie is an orthogonal 
projection map, we have 

=% (3.22) 

where denotes the closure of Y^ in the L^-topology. 

The initial data in Y^ has considerable gauge freedom, as only the gauge conditions 
Eq. 3.15 at B have been imposed, as well as (jP* = 0. Thus all smooth gauge transformations 
are allowed that preserve these conditions at B and approach asymptotic translations and/or 
rotations at inhnity. When we consider time evolution in Sec. 5, this gauge freedom would 
create a nuisance for dehning time evolution operators. Fortunately, it is possible to hx the 
gauge freedom completely as follows. The allowed gauge transformations SX := S*C*^ are 
precisely the ones lying in that additionally satisfy Eq. 3.14 and Eq. 3.15. Satisfaction 
of Eq. 3.14 requires that asymptotically approach a translation or rotation^^ at inhnity. 

Since M ^ 0 in the background, an arbitrary asymptotic Lorentz boost will yield 6Pi ^ 0. However, an 
arbitrary translation or rotation will not change M, Ja, or Pi to first order. Note that here the rotation 
may be arbitrary, i.e., it does not have to be along 
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From Eq. 3.16, it follows that satisfaction of Eq. 3.15 requires 


(3.23a) 

(3.23b) 

(3.23c) 


0 = & = 

0 = Modd = [31^3)A - 27r)!^0A^A + T^rT^rk “ ^ 

0 = (5h'even = “ {3^3A “ 27l^(j)f3A + '^r'^rk “ 

where M is the Ricci scalar of B. The requirements that ^-f^odd = Steven = 0 can be shown 
to imply ^|b = Crls = 0 by the same argument as given in [13] (with their f5^ = 2\/hn^ = 
27i^(f)^ and we have used the axisymmetry of the background only), which extended argu¬ 
ments given in [22] and [23]. Thus, the space Wg of allowed gauge transformations is given 
by 

Wg ■.={S*C*^ e 3*^ I ■CIb = ■Crlfi = 3a^'^\b = 0 and asymptotically 

® (3.24) 

approaches an arbitrary translation plus rotation at inhnity} 

Remarkably, Wg differs from only in that the condition = 0 in the dehnition of 
Wc has been replaced by the conditions ^\b = ^rls = 3 a^^\b = 0, and the asymptotic 
conditions at inhnity are somewhat different. 

We now hx the gauge completely on 3c by requiring orthogonality to in the L^-inner 
product. For smooth elements 6X G this L^-orthogonality will hold if and only if 


0 ;= CSSX = 0 (3.25a) 

((5M, 6P, (5j) (5(5X) = 0 (3.25b) 

(^5e,5zuAB'^ (‘55X)|b = 0 (3.25c) 

where 

SzUab ■= 3s^a {PB]r + T^r(l)kQB]c) , (3.26) 

and “(5J” in Eq. 3.25b includes all of the angular momenta, not just the ones conjugate to 
(j)\. The conditions at B expressed in Eq. 3.25c take the explicit form 

0 = fe(5(5X)|B = (3.27a) 

0 = 6zUABiS6X)\B = 3\^a {qB]r “ T^r ^t^^kP B]c) (3.27b) 


Thus, our gauge conditions on 5X correspond to the conditions Eq. 3.11, Eq. 3.14, and 
Eq. 3.16 with pab —t qab and gat —t —Pab except that the conditions requiring the vanishing 
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of the expansions, (5'(9odd and (5'deven, in Eq. 3.16b and Eq. 3.16c are deleted and replaced by 
Eq. 3.27b, and the condition 5J\ = 0 is replaced by 5J = 0. 

Let Yg denote the subspace of that is orthogonal to Wg. Then elements of Yg weakly 
satisfy our gauge conditions Eq. 3.25. Let fl^ : —)■ denote the orthogonal projection 

operator onto Yg. It is shown in Appendix A that, as for Lie, we have 

r- := Yg n . (3.28) 

Furthermore, as for Lie, we have that fl^ : Y^ is continuous in the (Frechet) topology 

of 

Now, by construction, we have Wg C since the gauge transformations in Wg were 
chosen to be precisely the gauge transformations in that he in Y^ . Taking perp-spaces, 
we hnd that Yg D {Y^)^ = (%)"*■, where the last equality follows from Eq. 3.22. This implies 
that (/ — Lie) = ng(/ — Lie) = {I ~ ne)!!^, which, in turn, implies that Lie and 11^ commute 

EeLlg = flglle. (3.29) 

The space of interest for us is 

Y^ := Ilglic[^^,] = Y^ n Y;^ . (3.30) 

This space contains a unique gauge representative of every element of Y^ . In the remainder 
of this paper, we will analyze dynamical stability for initial data in the space 


4. KINETIC AND POTENTIAL ENERGIES; POSITIVITY OF KINETIC EN¬ 
ERGY 

As discussed in the Introduction, we can use the (t-(;/))-reflection isometry, i, of the back¬ 
ground stationary black hole to decompose an arbitrary perturbation into its “odd” and 
“even” parts, P and Q, under the action of i. Let S be a (t-(;/))-reflection symmetric Cauchy 
surface for the exterior. We decompose the space of initial data, on S into L^-orthogonal 
parts Y°° = Y^^ © as follows. If the background is static, then i is purely a “t- 
reflection,” and we dehne P := (pafe,0) G Y^^ and Q := (0, g^b) G ^he stationary 

and axisymmetric case, we hrst decompose Pab and Qab into their “axial” and “polar” parts 
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with respect to the axial Killing fields 0^ as follows^^: 


Pab “1“ f3ab 'J ~ 0aA06© 

(4.1) 

Qab = ‘2af^(l)b)A + Pab + 

with 0^0^ = 0 = Pab = P{ab), Pab = P{abP PabpA = 0 = PabPX] and 7^® = 7^^®), 

= z/(^®l. Then the (t-0)-refiection odd and even parts of the perturbation respectively 

are: 

P '■= {Pab + 7^®0aA0b©, — ‘^0'^a4>b)A) ^ 'P^d 
Q ■= (2A()j0fe)A, Pab + V^^PaApbe) ^ 'P^en 

The negative sign in front of the a in the definition of P was chosen so that the transformation 
(p, q) I-© (P, Q) is canonical with respect to S defined in Eq. 3.8. 

Since the linearized constraint equations as well as our boundary conditions and gauge 
conditions are invariant under i, they cannot couple P and Q, so P and Q may be viewed 
as independent perturbations, i.e., if {P,Q) is a perturbation satisfying the constraints, 
boundary conditions, and gauge conditions, then (P, 0) and (0, Q) also satisfy the constraints, 
boundary conditions, and gauge conditions. 

The canonical energy is the map S’ : x —)■ R defined by 

S{5X, 5X) = ffs £t5X^ (4.3) 

Although the definition of S is not manifestly symmetric in 5X and 5X, it is, in fact, easily 
seen to be symmetric (see Prop.2 of [13]). An explicit formula for S can be found in Eq.86 
of [13]. Now, S is constructed from the background spacetime, so it is invariant under the 
reflection isometry, i, in the sense that for any perturbations SX and 6X in we have 
S{i*SX,i*5X) = S(SX,SX). It follows that under the decomposition P”" = Yp^^ © Ypp^, 
S cannot contain any (P-(5)-cross-terms. Thus, S splits up into two quadratic forms XP : 
'^odd X '^odd ^ K and ^ : Y^^n x ^ ® such that 

^[(P, Q), (P, g)] = ^(P, P) + ^(g, g), (4.4) 

where 

^(P,P) = ^[(P,0),(P,0)] (4.5) 

Note that this is a purely local decomposition into the parts that are parallel and orthogonal to the axial 
Killing fields. 
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^(4g) = ^[(o,g),(o,g)]. (4.6) 

We refer to ^ and respectively, as the kinetic energy and potential energy of the 
pertnrbation. 

In the simple case of a static backgronnd, where P = {pab, 0), it follows immediately from 
Eq.86 of [13] with pab = 0 that 

X = 2 (4.7) 

This expression is not manifestly positive dehnite since the second term in the integrand 
may dominate the first. However, pab is not arbitrary but is subject to the constraint 

D^Pab = 0 (4.8) 

as well as the boundary condition prr = 0 (see Eq. 3.16b). 

In the general case of a stationary-axisymmetric black hole, the kinetic energy is obtained 
by evaluating Eq.86 of [13] for a perturbation of the form 

Pab = I3ab + ; Qab = 2af^(t)b)K (4.9) 

Using the background constraints and integrating by parts using Eq. 3.16, we can write the 
kinetic energy in the form 

kXr = j N ]^{D^qabf - D^qabD^q’^^ +2 j N pl^- (p + ^ 

+ 'k k ^ + 2 ^ (p^^ + ^ Tjvgab • 

(4.10) 

The linearized constraints take the form 

Ca = D\b - 71^^ {2D[aab]A + a^Db^As) = 0 , (4.11) 

where Eq. 2.9 was used and, again, the form Eq. 4.9 for pab is understood. The linearized 
Hamiltonian constraint c = 0 is identically satished for perturbations of the form Eq. 4.9. 

We now rewrite Eq. 4.10 in a more explicit and useful form as follows. With the substi¬ 
tution Eq. 4.9, the hrst term in the integrand of Eq. 4.10 becomes 

je,= f N ^{D.qabY - D.qabD-q^’^ 

»/ s _ 

(4.12) 


23 



The last term of Eq. 4.12 can be written as 

2 [ = 2 [(a“-Zl'4>A=)p6«a0) + 2(a““D'<hAs)P[aa6]0)] 

Js Js 

(4.13) 

This expression can be simplihed by using the relation 

i?af,0A = (4.14) 

(which holds by virtue of (j)\ being Killing helds) and then eliminating Rab using the back¬ 
ground ADM equation Eq. 2.10 to obtain 

jja = -ANttSb ■ (4-15) 

Using this relation, we simplify Eq. 4.13 as follows: 

2 [ iV4>^®(a“=D'<hAH)(/^f,«a0) =-2 f [D'(iV4>^®D,$AH) «““«a0 + ^V=4^,<hA=] 

= 8 f iV(a“V) (7raA7r^)-2 f iV<l>^®(a“=D'$A=)(D,aae) - 2 f a“V0/^f,4>®“/^'<hA= 

V S V S V s 

= 4 [ N (tt.att^) - [ a“^aa04^f,4>®“/^'<hAH 

Je Je 

(4.16) 

Thus, we obtain 

je,=2 [ iV4>^® [(D[„a,]A + «fa4^b]4>A=)(4^f“«2 + a^'“4^''4>0r)l + 4 / iV (a“V) (^aA^r^) 

= 2 [ iV4>A0(D[,aJ)(D[Vl®)+4 [ iV (a“V) (^TaATT^) 

Je 7e 

To simplify the remaining terms in Eq. 4.10 we use 

'/hqca'^t = {o^a'^bA + a^TT^'®^^A^fe©) • (4.17) 

In addition, 

£Nqab = N^£^i^qab + 2qc(a(l>ADb)N^ = -4:Naf^7ib)A (4.18) 

where £^j^qab = 0 by axisymmetry, and Eq. 2.11 was used. We obtain 

££2 = 2 f N plf,- {P+^ qabT^""^^ 

+ 7^ ?c(a7rfe)) q^^'^T^d + I qcdKT^tq''^ + 2 ^ ^ 2g“'^7r^) 

= 2 ^ AT + (/© + 2af 7r®)“) ( 7 A© + 2a„(A7r^)) - ^ (/? + 7a + 2aM)^ 
-4^A^(a“^a®) (vTaATr©) 

(4.19) 
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Putting these results together, we obtain our final expression for 


^ = 2 ^ iV + ( 7 ^® + 2ai^7r®)“) ( 7 Ae + 2a„(A7r^)) 

(/^ + 7a+ 2aM)^ 


(4.20) 


We are now in a position to state and prove our first theorem: 


Theorem 1 (Positivity of kinetic energy). For axisymmetric perturbations of a stationary- 
axisymmetric black hole background, the kinetic energy Eq. 4.20, is a positive-definite 
symmetric bilinear form on 

Proof. It is instructive to first give the proof for the case of a static background, where i 
is given by a t-refiection and the kinetic energy is given by the much simpler expression 
Eq. 4.7. In that case, let f be the solution to the following boundary value problem (see [24] 
and [25]): 

= ; e|B = 0 ; e~0(l/p"-3)U . (4.21) 

Define P = {pab, 0) G by 

P := P -S*C* (4.22) 

where C* was defined by Eq. 3.12, i.e., we define 


Pab ■= Pab - DaDf,^ + hab^f, + Rabf, (4.23) 

where we have used the fact that for the static background, we have 7 r“^ = 0 and i? = 0 . 
Note that pab satisfies the constraint D'^pat = 0 and the boundary condition prrls = 0 (see 
Eq. 3.16b). Note also that p = = 0. We use Eq. 4.23 to eliminate Pab in favor of pab and f 

in Eq. 4.7. Integrating by parts and using the boundary conditions on (from Eq. 4.21) and 
prpB = 0, we find that the terms and the f-p cross-terms vanish. (This is essentially the 
same calculation as done in Eq. 3.18 and is an expression of the gauge invariance of JP with 
respect to gauge transformations that respect the boundary conditions Eq. 3.16.) Thus, we 
obtain 

JT = 2 y > 0 (4.24) 

Thus, JPf is non-negative and vanishes if and only if pab = 0. However, if pab = 0, then 
P = S*C* if, 0) e Wg. But P e C r C iT/. Thus, JT = 0 if and only if P = 0. 
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The proof proceeds in an exactly parallel manner for a general stationary-axisymmetric 

background, where is given by Eq. 4.20. We now let ^ be the solution to the boundary 

value problem 

(A - 27r^7r^) ^ ^ (/3 + = 0 ; ^ ~ O(l/p‘^-^)|oo (4.25) 

We again dehne P G by 

P := P-S*C* 0) (4.26) 

and we decompose P as 

^ = (/3ab + 7^®0aA0fee, -2af^(j)b)A'^ . (4.27) 

By a direct computation using the above dehnitions, we hud /3 + j^ + 2d^n‘l = 0. In parallel 
with the arguments in the static case we obtain (P, P) = (P, P). Hence, we obtain 

je = 2 j N [$Ae(/^[adJ)(P'“d'l®) + i3l, + ( 7 ^® + 2df 7 r®)“) ( 7 A 0 + 2d„(A7r^))] > 0 

^ (4.28) 

which shows that is non-negative. Furthermore, vanishes if and only if d^ = DaC^, 
Pab = 0 and 7 ^® = —2DaP^7r^^°'. This means ^ vanishes iff P = S*C* (O, he. 

P = S*C* {P -^^0a) e Wg. But P G PJd C r C So, = 0 if and only if P = 0, 
i.e., the kinetic energy is a positive-dehnite, symmetric, quadratic form on as we 
desired to show. 

□ 

Remarks: 

1. As shown in [16], a reflection isometry, i, exists for any subgroup of stationary- 
axisymmetric isometries that acts trivially and is such that its generators include 
the stationary Killing held and the horizon Killing held. If more than one choice of 
i exists, then Thm. 1 implies correspondingly stronger results. In particular, if one 
has a static, axisymmetric black hole, one has both a t-rehection and a (t-(;/))-rehection 
isometry. If we dehne and with respect to the t-rehection isometry, then the¬ 
orem 1 tells us that is positive dehnite. However, if we apply theorem 1 to the 
(t-0)-rehection isometry, we learn that, in addition, the potential energy, is also 
positive dehnite on “axial” metric perturbations, i.e., metric perturbations of the form 
Qab 20i^^(f)b)A- 
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2. The transformation P ^ P = P — S*C* (^, 0) is just a gauge transformation corre¬ 

sponding to making a normal displacement of the hypersurface S by The condition 
/^ + 7a + = 0 is just the condition that Stt = \/hp + qabT^°'^ = 0. Thus, writing 

the expression for ^ in terms of P corresponds to working in a gauge^^ where S is 
a maximal slice in the perturbed spacetime. The fact that JP{P, P) = {P, P) is a 

manifestation of the gauge invariance of 

3. Local Penrose Inequality: As shown in [13], positivity of the canonical energy S’ 
on is equivalent to the satisfaction of a local Penrose inequality. Consider pertur¬ 
bations of a Schwarzschild black hole with d < 7, with i chosen to be the t-reversal 
isometry. The fact that the Riemannian Penrose inequality holds [26] implies that 

is positive. However, since there do not exist nontrivial stationary perturbations of 
Schwarzschild with SM = SJ = 0 [27], [12], it follows that ^ is non-degenerate [13]. 
Thus, ^ is positive dehnite. Since we have just shown that is positive dehnite, it 
follows that S is positive dehnite. This implies that the Schwarzschild black hole is a 
strict local minimum of mass at hxed area of the apparent horizon. 

5. TIME EVOLUTION AND THE INVERSE KINETIC ENERGY HILBERT 
SPACE JP 

In this section, we will introduce the time evolution operators, /C and U, which will be 
seen to correspond to the quadratic forms, and ^, dehned in the previous section. We 
will then use the positive dehniteness of JS' to dehne a new Hilbert space that we will 
use in our stability analysis. The Hilbert space JS will not contain all of but it will 

contain all solutions of the form £tSX with 6X G 

The time evolution of the background with respect to P in Eq. 2.4 is a gauge transfor¬ 
mation generated by the lapse and the shift N_ = {N, A“), i.e. Eq. 2.4 takes the form 

X = S*C*N. (5.1) 

The evolution equations for smooth perturbations are obtained by linearizing this relation, 
i.e., 

SX = S*C*n + S*S{C*)N (5.2) 

Note, however, that this gauge choice is not compatible with the gauge conditions that we imposed in 
Eq.3.25. 
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where 


(5.4) 


n:^(SN,m‘‘) (5.3) 

and (5(£*) denotes the linearization of the operator C* given by Eq. 3.12. By a direct com¬ 
putation of S{C*) from Eq. 3.12, we obtain 

Qab ^nhab “1“ 2iV (j)ab d—lP ^o,b') ij^ab d—1^hab'j 

+ -^N (2gc(a7rfe) - ^qcdT^’"^hab - d^T^Qab) + £Nqab 

Pab 'j^achbd £ 2 ^ i^^qab 3“ Da^bq') 2 *^ ^ {Daqbc 3“ Dbqac Dcqab £ h^bDcq') 

- N {Rc(ab)dq'''^ + R\aqb)c - habR'^qcd) + habD'^ND^q^d - ND(^aD''qb)c 

- {NR - AN) gab + iNhab {-R^\d - Ag + D^D\d) 

+ iNq (7rac< - ^TrTTafe) - |iV [qcdK'^t + 2gc(a7rfe)d7r^‘^ - ^ {qcdT^''‘^T^ab + 2qc{aK)'^)] 

- ^ A [2pc(a<) - ^ {pT^ab + T^Pab)] “ l^q [Dc {N^TTab) “ 27rc{aD''Nf,)] 

Pc(a + '^qdc^{a A '^'^cqd{a D Nj,) — '^£Nqd{a'^b) 

(5.5) 


Dr 


N \Pab + ^Qd{a^b) 


-2 


+ [D^cNa)4 + D^rNb)7rt] 


where 


—A 7r“^ 
^nhab 


;= S*C* 



(5.6) 


is just the gauge transformation generated by n. In the above equations, we have used 
the background constraints and evolution equations, as well as the stationarity (but not 
axisymmetry) of the background. 

The perturbed lapse and shift n are to be chosen so that the solution lies in The 
linearized constraints Eq. 3.11 are preserved under time evolution, and thus C5X = 0 holds 
as an identity. Similarly, 5X identically satishes the boundary conditions Eq. 3.16 as can 
be checked by an explicit computation. Thus, 5X G and we only need to impose the 
conditions Eq. 3.25 so that 5X G . This requires that n satishes the equation 

CCn= -C5{C)K (5.7) 


This equation takes the form 

{d — l)A^n = ... 
2D^D(anb) = ... 
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(5.8a) 

(5.8b) 



where the dots denote terms of lower derivative order in n and source terms depending on 
5X. Thus, Eq. 5.7 is of elliptic character. The boundary conditions needed for the right 
hand side of Eq. 5.4-Eq. 5.5 to he in '7^°° are that S*C*n G Wg and that, in addition, n 
satisfy the boundary conditions arising from Eq. 3.25b and Eq. 3.25c. The conditions arising 
from Eq. 3.25c take the form 

0 = &(5(5X)|b = \{d - l)dln + l[{d- 2)^^ + lM]n + ... + 6e{6{C*)N) (5.9a) 

0 = 6wab{S6'X)\b = ^[AdrUB] + ... + dwAB{d{C*)X} (5.9b) 

where we have again suppressed the lower derivative terms in n and kept the source terms 
depending on SX in symbolic form. 

The problem of solving the above system for n corresponds precisely to applying the 
projection map^® Ilg to S*6{C*)T^. It follows from the results of Appendix A that a unique 
solution for n exists satisfying Eq. 5.7 and the above boundary conditions. Consequently, 
the evolution equations Eq. 5.4 and Eq. 5.5 take the form 

6X = S{6X). (5.10) 

where S : "7^°° —)■ is a linear map. Note, however, that since we have to solve the elliptic 

system Eq. 5.7 for n to determine the action of S on 6X, 8 is not a local differential operator. 
We have the following proposition: 

Proposition 5.1. Given any 5Xq G "7^°° there exists a unique solution SX{t) G to 
Eq. 5.10 that varies smoothly with t in the topology of and is such that (5X(0) = dX^. 

Proof. Uniqueness is obvious because the difference between two such solutions yields a 
solution to the linearized Einstein equation with vanishing initial data. Such a solution must 
be pure gauge, but the gauge has been completely fixed in so the difference between two 
such solutions must vanish. To show existence, we solve the linearized Einstein equation in 
linearized harmonic gauge to obtain a metric perturbation with initial data SXq. This 
solution will induce initial data SX{t) on each slice St. On account of the wave equation 
character of the linearized Einstein equation in this gauge, standard energy estimates (see, 
e.g.. Lemma 7.4.6 of [28]) can be used to show that SX{t) will he in and will vary 

It can be verified that S*S{£*)X lies in (see the proof of Prop. 5.1 below), so no He projection is 
needed. 
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smoothly with t. However, 6X{t) must also satisfy the linearized constraints Eq. 3.11 and 
the conditions Eq. 3.14 because SXq satishes these conditions and they are automatically 
preserved under time evolution. In addition “time evolution” (i.e., the evaluation of the 
solution on S* using the identihcation of with E given by orbits of maps B into B (see 
Fig. 1) and merely rotates it via the action of axial Killing helds. Since the perturbation is 
axisymmetric, it follows that Eq. 3.15 holds at all times if it holds initially. Thus, SX{t) G % 
for all t. The desired solution then is SX{t) = 'n.g5X{t). Smoothness of SX{t) in t follows 
from the smoothness of SX{t) in t together with the continuity of the gauge projection map 
Hg proven in Appendix A. □ 

The canonical energy, (S’, is related to the evolution operator S by 

(S{^, 6X) = Hs (<hA, S{6X)'^ = (Jx, SS{5X)^ (5.11) 

for all 6X,SX G This equation follows immediately from the dehnition of <S, Eq. 4.3, 
together with Eq. 5.10 and Eq. 3.6-Eq. 3.8. The relation Eq. 5.11 between cS and the time 
evolution operator S expresses the fact that (S is a Hamiltonian for the linearized equations. 

We now make use of the reflection isometry i of the background spacetime. As in the 
previous section, we decompose initial data (pab, Qab) G into its odd part, P G and 
even part, Q G under the action of i (see Eq. 4.2). Since the time evolution operator S 
is invariant under i, the gauge-hxed ADM evolution equations Eq. 5.4-Eq. 5.5 take the form: 


Q = K,P (5.12a) 

P = -UQ (5.12b) 


Here, the maps /C and lA act as 


/c: rjd ^ 


CXD 

even 


U : X 


OO 

even 


-yoo 
^ ^odd 


(5.13) 


We refer to /C and U, respectively, as the kinetic and potential time evolution operators. As 
in the case of 8, the operators /C and U are not local differential operators since the gauge 
hxing procedure involves solving an elliptic equation. 

We now return to Eq. 5.11 and decompose the perturbations SX and SX into their odd 
and even parts under i. Using the fact that no “cross-terms” can arise on account of the fact 
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that i is an isometry of the background spacetime, we see that the kinetic, potential and 
canonical energy quadratic forms, and are given in terms of the kinetic and potential 
time evolution operators, /C and U, by 


^(P,P) = (^P,SfCP^ 

^{Q,Q) = -{q,SUqY ( 5 - 14 ) 


Note that it follows immediately from the symmetry of the quadratic forms Jif and that 
both SfC : ^'^‘4 SU : —)■ symmetric linear maps in the L^-inner 

product. 

Taking the time derivative of Eq. 5.12a and using Eq. 5.12b, we obtain 

Q = -AQ (5.15) 


where 


A-.= KU-.%Zn^%Z.- (5.16) 

We wish to use spectral methods to solve this equation, but in order to do so, we need to 
dehne an inner product that makes A a symmetric operator. This can be achieved as follows. 
Let C denote the range of the operator fC. By Thm. 1, SfC is a positive-dehnite 

operator on 1^'^, so fC has vanishing kernel. Thus for all Q G /C[T(,‘^], there exists a unique 
P G that, Q = fCP. Using this fact, dehne a new inner product, (,)^ on /era 

by: 

(Q,g) ■.= J(r{P,P) (5.17) 

where P and P are such that Q = fCP and Q = fCP. That this is indeed an inner product 
follows from the symmetry, bilinearity and positive-dehniteness of (Thm. 1). Note that 
we can write this inner product in terms of the inner product Eq. 3.6 and the operator 
fC as 

(^, Q)^ = = (a ^Q) = p) (5.18) 

Thus, formally, the new inner product Eq. 5.17 corresponds to taking the matrix element of 
(iS/C)“^ in the inner product Eq. 3.6. 

We now complete the space fC[V^^ in the inner product (,)^ to obtain a Hilbert space 
Note that ■ However, does contain all Q G 11^® form 
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Q = JCP for P G In view of Eq. 5.12, this means that the even part of all perturbations 
that are of the form £t5X' for some perturbation 5X' G 1^°° will be represented in ^. As 
a consequence, our stability analysis of the next section will not apply to all perturbations 
in 1^°°, but it will apply to all perturbations of the form £t6X' with SX' G 

Finally, it will be convenient to complexify in order to apply spectral methods. We 
continue denote the complexified Hilbert space as i.e., we will not distinguish the com- 
plexihcation in our notation. 

6. DYNAMICS IN NEGATIVE ENERGY AND EXPONENTIAL GROWTH 

The operator A of Eq. 5.16 is well defined as an operator A : Jif —)■ with dense domain 

given by domA = For Q,Q E domA, we have 

{q,AQ^ {q,X{UQ)^ (sQMQ) = -{q,SUQ^ =^{Q,Q) ( 6 . 1 ) 

where Eq. 5.14 and the last equality of Eq. 5.18 were used. Thus A is densely defined 
symmetric operator on M'. Indeed, the entire purpose of introducing the inner product 
Eq. 5.17 was to make A symmetric. 

If AP were hnite-dimensional, then the operator A would be self-adjoint on and would 
admit an orthonormal basis of eigenvectors. This basis would diagonalize Eq. 5.15, and we 
could immediately conclude that any eigenvector of A with a negative eigenvalue would 
correspond to an exponentially growing solution. However, ^ is inhnite-dimensional and A 
is known only to be a symmetric operator defined on a dense domain in so the argument 
for exponential growth cannot be made so straightforwardly. 

Clearly A : ^ —)■ ^ is a real operator, i.e., it commutes with the (anti-linear) complex 
conjugation operator on ^ that we have just introduced via our complexihcation at the 
end of the previous section. We believe it likely that A satisfies the properties that it 
is bounded below and that it is essentially self-adjoint on the domain /C[lC“]. If these 
properties hold, then A would have a unique self-adjoint extension A that also would be 
bounded below. Unfortunately, it does not appear to be straightforward to establish either 
of these properties. Fortunately, as we shall see, we can bypass the issue of essential self¬ 
adjointness of A by simply choosing an arbitrary self-adjoint extension A] existence of such 
an extension is guaranteed by virtue of A being real. As we also shall see, we are able to 
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bypass the issue of A being bounded below by using spectral projections. 

Consider initial data of the form (Pq = 0, Qo = where Pq G Clearly, we have 

Qo G so Qo G Jf'. Let 5X{t) denote the unique solution in to Eq. 5.10 with 

this initial data (see Prop. 5.1). Let Qt G denote the (t-0)-even part of 5X{t). We 
claim first that Qt G ^ for all t. To show this, let 5X' be the unique solution in "P®" 
with initial data (Pq,(5o = 0). Then £tSX' has initial data (Pq = 0,(5o = ^Pq), and hence 
5X{t) = £ AX'it) for all t. It follows that Qt takes the form Qt = /CP/ where P/ G is 
the (t-(/))-odd part of 5X'{t). Thus, Qt G ^ for all t. Since dom^ = Cl JP, we also 
have Qt G dom^ for all t. Furthermore, when viewed as a one-parameter family of vectors 
in Qt satishes 

(see Eq. 5.15.) 

Let A be any self-adjoint extension of A on JP. Let {Pa} denote the family of projection 
operators associated to A such that^® 

OO 

A= j XdEx (6.3) 

— OO 

so that Ex projects onto the spectral subspace of ^ in (—oo. A]. Dehne Qt^jS by 


Qt,f} = (/ - E^)Qt (6.4) 

so that Qt^i 3 is the projection of Qt onto the spectral subspace {/3, oo) of A. Clearly, we have 
Qt = lim Qt^/s and \ \Qt,i3\\.^ < WQtWjt' for all A. Applying (/ - P^) to Eq. 6.2, we obtain 

/3—>—oo 


d^Qt,0 

dP 


—AQt,/3 ■ 


(6.5) 


Our next result will require the following lemma: 


Lemma 6.1. Let A be an arbitrary self-adjoint operator on a Hilbert space JP. Let Rt 
be a one parameter family of vectors in such that (i) Rt is twice differentiable in t, (ii) 
Rq = {dRt/dt)\o = 0 (iii) Rt lies in dom^l for all t, and (iv) we have 

^ = ( 6 . 6 ) 

for all t. Then P* = 0. 

The meaning of the integral (as a Lebesgue-Stieltjes integral) is explained, e.g., in Sec. 120 of [29] and in 
[30]. 
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Proof. Let A > 0 and let Rt^A be the projection of Rt onto the spectral interval [—A, A], i.e., 

Rt,A = {Ea — E-a) Rt (6-7) 

where Ex is the spectral family of projection operators associated with A. Then Rt^A satishes 
all of the properties of Rt stated in the hypothesis of this Lemma. Let 

/A(t) = II Rt,A f + II dRt^A/dt 11^ . (6.8) 

Then we have: 

^ = 23fJ {dRt,A/dt, Rt,A) + 23fJ {d‘^Rt,A/df, dRt,A/dt) (6.9) 

CL L 

= 2^ {dRt^A/dt, Ri^a) + 2‘St (—ARt^Ai dRt^A/d^ ( 6 . 10 ) 

where denotes the real part. Thus, we obtain 

^ II 7?t,A II + 2 II II + ARt,A (6.11) 

However, ARt^A < A^ || i?t A ||^ on account of the spectral projection Eq. 6.7. Thus, we 

obtain 

<(A2 + 2)/a (6.12) 

Since /a > 0 for all t G M and /A|t=o = 0, this inequality implies that fA(t) = 0 for all t, 

which implies that Rt^A = 0. However, the vanishing of Rt^A for all A implies that Rt = 0. □ 

Our instability results will be based upon the following proposition 

Proposition 6.1. Qt^/s as dehned by Eq. 6.4 is given by 

Qt,fi = cos(My^)(5o,/3 + hloQo,/? + cosh(Mi/^)(5o,/3 (6.13) 

where A+ = {I — Eq)A is the positive part of A, Hq is the projection onto the zero eigenvalue 
subspace^^ of A, and A- = — (Eq —Hoj^A is minus the negative part of A (so A^ is a positive 
operator). 

Proof. Let 

St,p = cos(My^)(5o,/3 + IfoQo,/? + cosh(Mi^^)(5o,/3 (6.14) 

IIo will be nontrivial if and only if there exist stationary perturbations for which 6M = 6J = 0. 


dfA 

dt 


dfA 

dt 
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Then is well defined since cos(t^y^) and Ho are bounded operators, and Qo,/? ^ 
donicosh(t^i'^^) on account of the spectral cutoff Eq. 6.4 used in defining Qt,/?- Since 
(5o,/3 £ dom^, it follows that St^p G dom^ for all t. Using the facts that dom^ C dom^Y 
and dom^ C doni^_ , it follows that that St^p is twice differentiable in t and 


d^Spp 

dt"^ 


-AS,^p. 


(6.15) 


Furthermore, we have S'o,^ = Qo ,/3 and {dSt,p/dt)\Q = 0. The proposition now follows imme¬ 
diately by applying Lemma. 6.1 to Rt = Qt,p — Spp- 

□ 


We now state and prove our main instability result: 

Proposition 6.2. Let Qq = /CPq with Pq G be such that ^{Qq,Qq) < 0. Then the 
solution generated by the initial data (Pq = 0, Qq) grows exponentially with time in the 
sense that there exists U > 0 and a > 0 such that 


WQtWj^ > C (6.16) 

Proof. By Eq. 6.1 we have 

{Qo,AQo)j^ = '^{Qo,Qo) <0. (6.17) 

Therefore, there exists a > 0 such that E_a 2 Qo ^ 0. Clearly, we can choose fd sufficiently 
large and negative so that P_q 2 (/ — Ep)Qq = E_a 2 QQ p ^ 0. However, by Prop. 6.1 and the 
spectral representation Eq. 6.3, we have 

II P-Q2Qt,/311^ = / cosh^(tA^/^)d|| PaQo lljr 

J-o? 

> cosh^(ta) II E_^ 2 Qq^p 11^ > exp(2at) || P_q2Qo,/ 3 Hjf’ (6-18) 

Thus, we obtain 


II Qt \\j^ ^ II Qt,P llj^ P II E_QpQt,p llj^ ^ II P—||j^GXp(Q;t) 

as we desired to show. 


(6.19) 

□ 


Remark; By similar arguments, it follows that 
However, we have 



also grows exponentially with t. 


Qt 


= II /CP 






= jr(R,Pt) 


( 6 . 20 ) 
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where Pt is the (t-0)-odd part at time t of the solution generated by the initial data (0, Qo). 
Thus, the kinetic energy, of the solution generated by (0, Qq) blows up exponentially 
with time. Since the total canonical energy S’ is conserved, it follows that —^{Qt,Qt) also 
blows up exponentially. Since JS and are gauge invariant, this shows that the exponential 
blow-up found in Prop. 6.2 is not a gauge artifact. 

We now reformulate our results as a Rayleigh-Ritz-type of variational principle for ana¬ 
lyzing black hole instability, which makes no direct reference to the Hilbert space 


Theorem 2 (Variational Principle for Instability). For any P G consider the quantity 


co^P) 


^{fCP,fCP) 
JP{P, P) 


( 6 . 21 ) 


Then if < 0, the solution SX{t) determined by the initial data (P, 0) will grow with time 
at least as fast as exp(Q;t) for any a < |a;|, in the sense that the kinetic energy, JP, of £tSX 
will satisfy 

lim [jP{£t6X, £t6X) exp{—2at)] = oo . ( 6 . 22 ) 

t^OO 


Proof. Let Qo = /CP. Then we have 


^(/CP,/CP) = ^(go,Qo) = {Qo,AQo)j^ 


(6.23) 


and 


Thus, 


^(p,p) = (go,go) 




2 _ (go, -^go) 

CO — 


(6.24) 


(6.25) 


(go, go) 

so if < 0 , then it must be the case that E_cfiQo 7 ^ 0 for any o£ < 

The solution £t5X has initial data (0, go). The theorem follows immediately by applying 
Prop. 6.2 and the remark below that proposition to this solution. 


□ 


Remark; The variational principle in Thm. 2 may be viewed as a generalization of the 
variational principle of [31] for spherically symmetric perturbations of static, spherically 
symmetric spacetimes (with matter helds), which itself is a generalization of the variational 
principle of [32]. Indeed, much more generally, in an arbitrary theory, the key ingredients 
needed to obtain such a variational principle for linearized perturbations off of a stationary 


36 



background are that (i) the theory be derived from a Lagrangian, so that the canonical 
energy, ^^ of perturbations is well dehned; (ii) gauge conditions can be chosen that nniqnely 
hx the gauge in snch a way that S acts as a Hamiltonian for the linearized perturbations and 
time evolution yields a snitably smooth map on the initial data space; (hi) the background 
admits a snitable reflection isometry, i, so that the canonical energy S can be decomposed 
into a kinetic part and a potential part ^; (iv) the kinetic energy is positive dehnite. 
Thus, it should be possible to generalize onr results to a wide variety of other cases as well 
as a wide variety of other theories. As simple examples, in Appendix B we analyze the 
cases of an axisymmetric Klein-Gordon scalar field and an axisymmetric electromagnetic 
held propagating in an arbitrary stationary-axisymmetric black hole backgronnd with a 
(t-0)-rehection isometry. 
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Appendix A: Projection Operators 

The purpose of this Appendix is to outline the proof of the following proposition: 

Proposition A.l. The orthogonal projection operators He : ^gr —)■ and H^ : —)■ 

dehned in Sec. 3 map to itself and—as linear maps on —are continnous in the 

natural (Frechet) topology on 

The basic strategy for proving this proposition consists of constrncting these projection 
operators by solving elliptic equations. The nice properties of the projection operators then 
follow from elliptic regularity. This strategy was successfully implemented in Proposition 5 
of [13], with the main technical steps carried ont in Lemma 3 of that reference. However, 
the choice of boundary conditions for spaces onr Wc and Wg diher from the space, WhW) 
considered in [13], and there are several other (relatively minor) diherences as well as some 
(small) improvements that can be made. Therefore, rather than merely referring the reader 
to [13], we shall we shall now ontline the main steps needed to prove Prop. A.l. However, 
our proof mirrors the proof given in [13] in its essential details. 
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To project SX G ^gr orthogonally to with dehned by Eq. 3.17, we would like 

to solve the elliptic equation 

CC*^ = CSX (A.l) 

with ^ satisfying the conditions appearing in the dehnition Eq. 3.17 of and also such that 

SX — C*^ satisfies the boundary conditions Eq. 3.14 and Eq. 3.15. If we can hud such a 

then 

SX' := SX - C*^ (A.2) 

will satisfy CSX' = 0 together with the boundary conditions Eq. 3.14 and Eq. 3.15. Thus, 
we have SX' G % whereas S*C*^ G Wc, so the desired projection map is 

U^SX = SX' = SX- C*^. (A.3) 


Similarly, to obtain the projection operator Ilg, we we wish to solve the equation 

CC*^ = CSSX (A.4) 

with ^ satisfying the conditions appearing in the dehnition Eq. 3.24 of Wg and also such that 
SX — S*L*^ satishes the boundary conditions Eq. 3.25b and Eq. 3.25c. We then obtain 

EgfjX = (A.5) 


Thus, the proof hinges on the ability to solve the equation CC*^ = a for suitable a such 
that ^ satishes the desired boundary conditions. The key step in showing this is the proof 
of a Poincare inequality of the form 


I i \\Wf(BWj, 


< CMpi 


Pk llL20i2 — ^ 


i IIW2®Wpi 


(A.6) 


for positive constants c, C where Wp®Wp denotes the weighted Sobolev space on .^ = (,^, 
and the weighting matrix appearing in the middle term is 


Mp-.= 



(A.7) 


where p is a positive function that goes to 1 in a neighborhood of B and approaches 
(xf + ... + near inhnity. We want this inequality to hold for all ^ lying in a suit¬ 

able closed subspace of © Wp. As can be seen from the proof of lemma 3 of [13], a 
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certain boundary term at B arises in the estimates, and, for the validity of the first inequal¬ 
ity in Eq. A. 6 , we need ^ to be such that this boundary term vanishes. This requires the 
vanishing of 

[ r“ {D.U’’ - - irDan 

J B J B 

= [ i^rdr^r + " ^rdr^r + (A.8) 

J B 

= f = "2 / Ai^ 

J B J B 

Thus, the Poincare inequality Eq. A .6 will hold if we restrict ^ to the closed subspace of 
IT^ © Wp dehned by^® {,^Ai^)\B = 0 (or alternatively, the closed subspace dehned by 

= 0). In addition, we may impose the vanishing at B of any other components of ^ 
and and/or their tangential derivatives. However, we cannot impose conditions on radial 
derivatives of these quantities, since such conditions would not dehne closed subspaces of 
W^p © W}. 

The closed subspaces of Wp © Wp of interest are 

^HW ■■= (i = (e,r) e Wl © Wl I eis = e.ls and = O} (A.9a) 

■■= (i e Wl © Wl I = 0} (A.9b) 

^g-={i^Wl®Wl I e|B=er-|B = («^)|B=0} (A.9c) 

The space ^ will be used for the construction of He and the space ^ will be used for the 
construction of H^. The space ^hw is included in this list because it is the space used in 
[13], but it will not be used here. 

As shown in the proof of lemma 3 of [13], for d > 4 (i.e., spacetime dimension 5 or 
greater), the Poincare inequality Eq. A .6 holds on the spaces ^hwi and ^g. Equation 
Eq. A .6 also holds in d = 3 provided that, if the black hole is non-rotating (i.e., for the 
Schwarzschild case), we pass to the subspace orthogonal in to where is 

the timelike Killing held of the background black hole. 

Now, let denote any of the spaces ^hw, and with the additional projection 
orthogonal to taken in the case of a Schwarzschild background in d = 3. In view 

This condition is well defined since for any G W/, its restriction to B is well defined in L'^{B) and its 
tangential derivatives are therefore also well defined (weakly). 
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of Eq. A.6, the Lax-Milgram theorem^® implies that given any [3 G there exists t] & 
such that for all y G we have 

{CM,2,CM,x)^{^,x) , (A.IO) 

where (,) denotes the inner product. In particular, since this equation holds for all 
X G it follows that 

e := (A.ll) 

is a weak solution to 

CC*^ = a (A.12) 

where a = Aipl3. However, since CC* is elliptic (see [24], [33] and also Ch.VII Sec.2 of [21]), 
if ^ G ^ Hfc {W^ © Wj;) (so that a G (Mp^) {p‘^W^ © pWj^)), then ^ e {Mp^) 

{p^Wp © pWp) and Eq. A.12 holds in the ordinary (strong) sense. Furthermore ^ depends 
continuously on a in the natural (Frechet) topology of these spaces. It also follows from 
Eq. A.IO together with the fact that Eq. A.12 holds in the strong sense that for all 
(j\4p^) Dk (p^Wp © pWp), the solution ^ to Eq. A.12 given by the Lax-Milgram theorem 
satishes 

{C(,Cl) = {CC(,l) (A.13) 

However, since ^ and ip are smooth, the difference between the left and right sides can be 
explicitly computed by the same calculation as previously done in Eq. 3.18. In this case, no 
boundary terms arise from inhnity and we obtain 

{C*^, C*^) - {CC*^, ip) = ‘^ [ {-dri^ 6e +Ip (5^9even + di3odd “ Szuab) (A. 14) 

J B 

where Szuab is given by Eq. 3.26 and we have used ^Ai’^ = 0 to replace ip"^ by for 

some ip^^ = ip^^^i Since the right side vanishes for all ± G {Mp3^) Hk {p^W^p © pW^p), we 
obtain the following additional boundary conditions on the solutions ^ in the various cases: 

• 3^hw- Definition: = 0; Additional conditions on solutions: 5e\B = 0, 

SiP+Ib = 0 . 

• Definition: = 0; Additional conditions on solutions: Se\B = 0, ^devenls = 0, 

^'doddlB = 0. 

This theorem is just the Riesz lemma applied to (/3, •), which is a bounded linear map on Wp, © Wp and 
hence, by Eq. A.6, is also a bounded linear map on the inner product space defined by {C*Mp •, C*Mp ■). 
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• ^g\ Definition: = 0 = = 0; Additional conditions on solutions: 

dels = 0, 5tuab\b = 0. 

Thus, the solutions given by the Lax-Milgram theorem to Eq. A.12 with a G Dk 

(p^Wg © pWg) satisfy precisely the conditions at B needed to obtain the spaces Whwi 
and Wg, respectively. 

The key ingredients are now in place to construct the projection map lie. As explained 
at the beginning of this section, to construct lie, we wish to solve Eq. A.l for any given 
5X G In order to obtain a solution that satisfies the desired boundary conditions at 

B, we choose ^ be smooth and of compact snpport such that: (i) ^^\b = 0; (ii) C*^ — SX 
satishes the horizon bonndary conditions Eq. 3.15 reqnired for elements of Xp, and (hi) 
— C5X^ |_B = 0. There is no difficulty simultaneously satisfying all of these condi¬ 
tions because ^ is not constrained by any eqnations holding in the bnlk, so all of the various 
radial derivatives at B can be chosen independently. 

We now solve 

CC*l = C5X - CC*i (A.15) 

Condition (hi) gnarantees that the right side lies in pWp)i so we obtain a 

solution G {M.pt^c)t^k{,P^Wp®pWp). It follows that (5X — £*(,f+,^) satisfies the constraints 
Eq. 3.11 and the horizon bonndary conditions Eq. 3.15 needed to be an element of 
However, it need not satisfy the remaining conditions Eq. 3.14, namely 5M = SJa = dPi = 0. 
Since these are only a hnite nnmber of conditions, only a hnite co-rank further projection is 
needed. 

In [13], this additional projection was done by choosing an arbitrary minimal set of 
solutions T/ that satisfy the horizon bonndary conditions and whose span yields arbitrary 
values of {dM,6J\,dPi). The hnal projection to was done by subtracting from 6X — 
+ the linear combination of T/ that makes dM = 6J\ = dPi = 0. However, this does 
not, in general, yield an orthogonal projection, as would be needed to obtain commutativity 
with Hg. 

We will carry out this additional projection by making a special choice of T/ that yields 
an orthogonal projection. Let be a basis for asymptotic translations and rotations with 
respect to the axial Killing helds at inhnity, with each chosen so that it vanishes in a 
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neighborhood of B. Now solve the equation 

CCXj = -CCx, (A.16) 

to obtain a solution in (MpS’c) Dk {p^W^ © pW^). Dehne 

ij = Xj + Xr (A. 17) 

It follows immediately from the dehnition Eq. 3.17 of that S*C*^^ G #);. Any element 
of must therefore be L^-orthogonal to C*^^. In particular, we have ^ But 
satishes the constraints Eq. 3.11 and the horizon boundary conditions Eq. 3.15 needed to 
be an element of '7^°°. This implies that cannot satisfy 6M = 6Jx = SPi = 0. More 
generally, no linear combination of the can satisfy SM = SJ^ = 5Pi = 0, so the span 
of must yield all possible values of ((5M, (5Ja, 5Pi). Thus, may be used to play the 
role of T/. 

Now let 

i = l + i+E^4 / (A-18) 

I 

where the c/ are chosen so that 5X — satishes 5M = 5J\ = 5Pi = 0. Putting everything 
together, we have shown that any 5X G can be written as 

5X = 5X' + C*i (A.19) 

where S*C*^ G and 6X' G Thus, we have IIc^X = SX' G as we desired to 

show. The continuous dependence of SX' on SX in the topology of follows directly 
from the continuous dependence of the solution on a in Eq. A.12. 

The corresponding results for the projection map 11^ are obtained in complete parallel 
by replacing and with Wg and ^ and inserting the map S in appropriate places. 

Appendix B: Scalar and Electromagnetic Fields 

In this Appendix, we give a treatment of scalar helds and electromagnetic helds on a black 
hole background analogous to our treatment of gravitational perturbations. The black hole 
background is assumed to be stationary and axisymmetric with a {t-4>) rehection isometry i 
and with a bifurcate Killing horizon with bifurcation surface B, but it need not be a solution 
to Einstein’s equation. As in the gravitational case, we consider only axisymmetric scalar 
and electromagnetic helds. 
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1. Scalar Fields 


The case of a Klein-Gordon scalar field 0 propagating on a black hole background is 
particularly simple because there are no constraints and there is no gauge freedom. The 
initial data for a scalar held is SX = (vr,^), where 0 = 0|s and n = (m^V^0)|s, where is 
the unit normal to S. The analog of is 

^KG:=i^'(S)©T2(S) (B.l) 

with the natural volume measure on E. The analog, of is the intersection of the 

corresponding weighted Sobolev spaces (see Eq. 3.4). The symplectic form is 

:= ^(^0-07r) , (B.2) 

which is represented on by the orthogonal linear map S given by 


S{7T,ii) = 


(B.3) 


Since we have no constraints or gauge conditions to enforce, the analogs of the projectors 
flc and Ilg are trivial, so 

In both the static and rotating cases, the rehection odd and even parts, respectively, of 
the perturbation are simply P = (tt, 0) and Q = (0,0). Note that in this case we have 
— ^^^even- The kiuetic and potential energies for axisymmetric helds are: 


— / N7r‘^ 
JT. 


"^KG = N [{Daijf + 


(B.4) 


Both JPkg and '^kg ar® manifestly positive dehnite. Indeed, for a scalar held^'^, the canonical 
energy takes the form 


f^KG — ’^YJG + ^¥JG — 2 / 


(B.5) 


where is the energy-momentum tensor 




(B,6) 


20 


For other fields the canonical energy can differ from the right side of Eq. B.5 by a “boundary term”; see 


the appendix of [34] for further discussion. 
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For axisymmetric fields we then have = NT^yU^u’^ > 0, which again shows the 

positivity of the canonical energy. 

The evolution equations for axisymmetric helds are 


ij) = Ntt 

71 = Da {ND^Tp) - iVmV 


(B.7) 


Thus, the operator /Ckg corresponds simply to multiplication by N. It is then obvious that 
for tjj G I^i^even) have tjj G A^Kcl^^odd] if if '’P\b = 0. The “inverse kinetic 

energy Hilbert space” is then the L^-space on S with volume measure given by N~^ 

times the natural volume measure on S. 

Since '^kg is positive dehnite, the operator ^kg : ^kg —t J^kg is positive dehnite, so 
there cannot exist any exponential growth instabilities for axisymmetric^^ scalar helds in 
any black hole background. In fact, in this case the expression for ^kg is sufficiently simple 
that one can prove stability for elements of l^^even ff^^f ii® i'^ '^^G; as done in [2] for the 
case of Schwarzschild. The restriction to data vanishing at B can then be eliminated (for 
general black holes) by the “trick” used in [5]. For the case of a Schwarzschild or Kerr black 
hole, these results can be greatly improved by the methods of [3] and [9]. 


2. Electromagnetic Fields 

The initial data for a Maxwell held is 5X = [Ea, Aa), where Ea is the electric held on S 
and Aa is the pullback of the vector potential to S. The analog of ^gr is 

^em :=f^"(S,T*S)©L2(S,T*S) (B.8) 

and the analog, of is the intersection of the corresponding weighted Sobolev 

spaces. The symplectic form is 

<5X) := [ {E^Aa - AaE ^), (B.9) 

./e 

which is represented on by the orthogonal linear map S given by 

S{E\Aa) = {A%-Ea). (B.IO) 

For a static black hole background, the restriction to axisymmetric perturbations is unnecessary. 
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The analog of the operator C of Eq. 3.10 maps smooth initial data on S to smooth 
functions on S and is given by 


C 


EM 


Ea 


;= D^E„ 


(B.ll) 


i.e., the constraints are 

Cem = = D°'Ea = 0. (B.12) 

The formal adjoint of £em maps smooth functions on S to smooth initial data on S and 
is given by 

^ -Dai 


'^EM'C — 


0 


(B.13) 


The gauge transformations are then 5^X = 
The analog of the space Wc of Eq. 3.17 is 


(^em)c := e I eii? = 0 and ^ ~ const.U} (B.14) 

The symplectically-orthogonal space is (Tem)c := If SX G then SX G 

{XEMjT ■= ('% m)c n if and only if 


Cem — 0 

Q = 0 


(B.15) 


where 


Q: = 


Ea 


(B.16) 


J OO 

is the electric charge of the solution. In view of Eq. B.12, the charge integral may be taken 
over any sphere homologous to a sphere at inhnity. 

The analog of the space Wg of Eq. 3.24 is 


e ^em} 


(B.17) 


i.e., ('pEM)g consists of all gauge transformations lying in The space (Tem)^ is the 

orthogonal complement of Elements of (Tem)“ := satisfy the gauge 

conditions 


0EM •— — D'^Aa — 0 

Arls = 0 . 


(B.18) 
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These conditions fix the gauge completely. The projection operators (nEM)c and (nEM)c/ can 
now be defined in direct analogy with He and 11^ and they satisfy the analogous properties. 
Indeed, in the Maxwell case, the operator /Iem-^em = is just the Laplacian operator 
acting on scalar functions, so the analogs of many of the results of Appendix A are standard. 
The perturbations in 

•= (B.19) 


can again be split into their odd and even parts under the reflection symmetry i. In the case 
of static black hole background, we may take i to be the t-refiection isometry. The reflection 
odd and even parts of the initial data are then P = (i?a, 0) and Q = (0,Aa), respectively. 
The kinetic and potential energies are given by 


P^EM — J 

^EM = 2 


NEi 




(B.20) 


Both and '^em are manifestly positive-definite for any static black hole background. 
Thus, as for scalar fields, there cannot exist any exponential growth instabilities for electro¬ 
magnetic fields propagating on an arbitrary static black hole background. 

However, the situation is quite different for axisymmetric electromagnetic fields propa¬ 
gating on a stationary-axisymmetric black hole background. We can decompose Ea and Aa 
into their 0-refiection odd and even parts as 


Ea — E (pXa + Ea 


(B.21) 


Aa — A^(j)\a + Aa 


(B.22) 


with E^(j)\a = A'^(f)\a = 0. The reflection-odd and even parts of the initial data are then, 
respectively, 

-PeM = (^Ea, —A^(j)\a^ (B.23) 

Qem= (E^(j)Aa,Aa) (B.24) 


The kinetic and potential energies are then given by 


— 




NEl + 2N 


= N 


E^E\ + 2(D^aAb])^ 


(B.25) 
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It is clear that there exist black hole backgrounds for which J^m fails to be positive^^. To 
see this, we simply start with any black hole background spacetime (e.g., Schwarzschild) and 
any initial data {Ea, Aa) G ^^^ch that (j)Q&D[a{A^ 7 ^ 0 at some point a; G S. Then 
we can choose a function iV'® that is sufficiently large in a neighborhood of x so that if we 
replace iV® by iV'® then the last term in the expression for will be negative and will 
dominate the first two terms. We can then construct a new stationary-axisymmetric black 
hole spacetime that has the same induced spatial metric on S, but with a new stationary 
Killing field that has the same lapse function N as but has shift vector On 

this new black hole spacetime, will be negative for the original initial data {Ea, Aa). Of 
course, the new stationary black hole spacetime will not be a solution to Einstein’s equation. 
It is not obvious whether can be made negative for any black hole background that is 
a solution to Einstein’s equation. 

The evolution equations for the Maxwell held can be written as: 


Aa — —Da(p + NEa — 2N^ D\^aAh] 

Ea = ‘2D’’ [ND[bAa] + N[bEa]) 

where the gauge conditions Eq. B.18 applied to 6X imply that ip satishes: 


(B.26) 


-Ap + D^{NEa - 2N’’D[aAb]) = 0 
{-drP - 2AN’’D[aAb])\B = 0 


(B.27) 


The operators JCem and Uem can be read off from these equations by taking their even and 
odd parts under the rehection symmetry i. 

Since need not be positive dehnite, we cannot directly apply the results of Sec. 5 
and Sec. 6 to the electromagnetic case. Remarkably, however, the potential energy '’^em is 
positive dehnite (see Eq. B.25). Therefore, if we simply make the canonical transformation 
(Tem, Qem) {Pem = Qem, Q'em = —Pem), then all of the results of Sec. 5 and Sec. 6 apply 
(-Pem)Q™)- particular, if SX G and £tSX has negative kinetic energy, then SX 
grows exponentially with time. 


This possibility was first pointed out to us by A. Ishibashi. 
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